Enhancing single-parameter quantum charge pumping in carbon-based devices 
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We present a theoretical study of quantum charge pumping with a single ac gate applied to 
graphene nanoribbons and carbon nanotubes operating with low resistance contacts. By combining 
Floquet theory with Green's function formalism, we show that the pumped current can be tuned 
and enhanced by up to two orders of magnitude by an appropriate choice of device length, gate 
voltage intensity and driving frequency and amplitude. These results offer a promising alternative 
for enhancing the pumped currents in these carbon-based devices. 
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A current flow in a two terminal conductor normally 
requires a bias voltage difference between the electrodes. 
However, in systems where the coherence length of the 
electrons exceeds the device length, the use of ac fields 
can produce a non-vanishing current even at zero bias 
voltage. This quantum coherent effect is called quantum 
charge pumping [J-Q and has been observed experimen- 
tally When the driving frequency is low enough 
in such a way that f2 <C (being r the time spent by 
an electron to traverse the sample) then the transport is 
said to be adiabatic. Under such conditions at least two 
out of phase time-dependent parameters arc necessary to 
achieve pumping Q . Beyond the realm of the adiabatic 
approximation, pumping with a single parameter is also 
possible (sl-fl^ as confirmed by recent experiments [Tsl - 
Il5| . While a drawback for such monoparametric pumps 
is that higher frequencies are required to achieve similar 
currents, a big advantage comes from the reduction in 
the number of necessary contacts which makes it promis- 
ing for scalable, low-dissipation, devices. In this way, 
single-parameter pumping has become a quest for achiev- 
ing higher currents in nanoscalc transmission channels. 

The advent of graphene [l6| and their lower dimen- 
sional cousins, graphene nanoribbons and carbon nan- 
otubes , provided an outstanding ground for exploring 
the physics of quantum charge pumping as evidenced by 
recent studies. Most of them are focused in flat graphene 
samples driven in the adiabatic regime [isl - Eoj or beyond 
the adiabatic regime but in the energy range close to the 
Dirac point [2ll.l22|. The observation of other prominent 
features such as van Hove singularities usually requires 
an important gate voltage but recent experiments per- 
formed on twisted graphene layers demonstrate that sin- 
gularities can be brought arbitrarily close to the Fermi 
energy f23j , opening interesting prospects for beyond-the- 
Dirac-point-studies. 

In this Letter, we explore the interplay between elec- 
tronic structure and non-adiabatic effects in carbon- 



based quantum pumping devices with a single parameter. 
The transmission channels are composed of carbon nan- 
otubes or graphene nanoribbons but our results are ex- 
pected to be valid for generic quasi one-dimensional sys- 
tems. By modeling a device with low resistance contacts 
and driven by a single time-dependent gate, we show that 
the pumped current can be enhanced by up to two or- 
ders of magnitude by gating the system close to a van 
Hove singularity. Moreover, it is shown that substan- 
tial improvement is also possible at the Dirac point when 
tuning the interplay between driving frequency and de- 
vice length. Our study points out alternative directions 
to bring these devices closer to reality. 

Tight-binding model and Floquet solution. We consider 
a device consisting of a graphene nanoribbon with passi- 
vated edges, or a single wall carbon nanotube, of length 
L connected to two semi-infinite electrodes. For the sake 
of simplicity, the electrodes are considered to be a pro- 
longation of the sample located at the centre. The sys- 
tem is described through a single 7r-orbital Hamiltonian 
M He = J2i Eicfc, - J2{i,j) [li,j4cj + H.C.], where c+ 
and Cj are the creation and annihilation operators for 
electrons at site i, Ei are the on site energies and 7ij- 
are nearest-neighbors carbon-carbon hoppings. To sim- 
ulate a system with low resistance contacts, the central 
part of length L (the "sample" ) is connected to the semi- 
infinite leads through matrix elements 74 smaller than 
the carbon-carbon hopping 7 = 2.7 eV We 
set 74 = 0.77 in order to mimic the Fabry-Pcrot con- 
ductance oscillations [25| experimentally observed at low 
temperatures [26| . 

In addition to the possibility of homogeneously gating 
the whole central sample, which is useful for revealing 
Fabry-Pcrot oscillations as in [2^, here we include an 
ac gate asymmetrically disposed on the scattering region 
as represented at the bottom of Fig. 1. The ac gate, 
which is responsible for the dynamical breaking of the 
left-right symmetry, is modelled through an additional 
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FIG. 1: (color online) Inset: Pumped current as a func- 
tion of the gate voltage for a (24, 0) carbon nanotube with 
L — 100 nm. Main frame: maximum current pumping as a 
function of the system length at two fixed driven frequencies: 
(solid red line) hQ,2 = 4.050 meV and (dotted line) ftfii = 
2.025 meV. The ac field intensity is set to eVac = 1.35 meV. 
The vertical down arrows (grey) mark the points where the 
pumped current vanishes. Bottom panels: schematic repre- 
sentation of the system with symmetric (left) and asymmetric 
(right) ac gate. 



on site energy term Ej^Q = eV^c cos{nt) + eVgate where 
G is the region where the ac gate is applied, Vac is the 
amplitude of the ac gate, and Vgate is the one of a static 
gate which is applied to the whole sample, fl is the driven 
frequency. The current is obtained by using Floquct the- 
ory [1, [23, [13 • The dc component of the time-dependent 
current I{t) can be computed as 
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where T^J_j^{s) are the transmission probabilities of the 
carriers coming from the left (L) to right (R) electrodes 
which might absorb or emit | n \ photons depending 
if n > or n < 0, respectively. These transmission 
probabilities can be fully written in terms of the Flo- 
que t Green's functions for the system as explained in 
The number of Floquet channels (n) used in the 
calculation is chosen by requiring the convergence of the 
transmission functions with at least six significant digits. 
The electron distribution in the left (right) lead is given 
by the Fermi function /i (e) (/^Xe)). The non-interacting 
model is justified when screening by the surrounding gate 
or a metallic substrate reduces electron-electron interac- 
tions, ruling out effects beyond our present scope [2^. At 
the same time, electron-phonon interactions can be ruled 
out at low bias voltages (< 150 meV) since the relevant 
inelastic mean free path for scattering with acoustical 



phonons is on the order of 1-2 micrometers at 300 K 130| . 
Recently, a similar method was applied to unveil laser- 
induced band gaps in graphene and their effect on charge 
transport 31 1. 

Single-parameter pumping: interplay between elec- 
tronic structure and non-adiabatic effects. Although in 
the following we present results for metallic nanotubes, 
the same physics was found for armchair edge graphene 
nanoribbons. This is consistent with the argument pro- 



posed in [32| based on the use of a mode decomposition 
scheme for the electronic states. Furthermore, we expect 
the trends found here to be valid for a generic quasi one- 
dimensional system. Figures 1 and 2 show the maximum 
pumped current calculated for a (24, 0) carbon nanotube 
with an asymmetrical ac gate (see bottom schemes in 
Fig. 1) extending over a typical length Lac = 30 nm. 

When pumping close to the charge neutrality point 
(CNF), the Fabry-Perot oscillations manifest as an oscil- 
latory pumped current as shown in the inset of Fig. 1. 
Its period is commensurate to the energy level spacing, 
A, which characterizes the system spectrum. Changing 
the device length allows the tuning of the energy level 
separation since A cx 1/L. On Fig. 1 (main frame), 
one can notice that the maximum pumped current can 
be strongly modulated by the device length. The solid 
red line corresponds to a driving frequency which doubles 
the one for the dotted line. The zeroes in these curves 
(marked with down arrows) can be understood in terms 
of a wagon-wheel effect [12, [11]: whenever the driving 
frequency is commensurate with the level spacing, the 
system behaves as in the static regime and no current is 
pumped through the system. A similar phenomenon was 
noted in [^ for a double barrier structure. Another fea- 
ture in the decrease of the successive maxima in Fig. 1 
(main frame). To explain this issue, first one notes that 
the pumped current is obtained by integrating a kernel 
function given in Eq. [TJ which in our case turns out to 
be periodic as demonstrated on Fig. 1 inset. While the 
kernel amplitude remains approximately constant when 
changing the device length, the relevant energy scale over 
which the integration is performed is min(hn, A). This 
leads, for big enough L to a series of maxima which scale 
with A oc 1/L as seen in Fig. 1 (main frame) . 

Beyond the CNF, application of a gate voltage can un- 
veil interesting features when tuned nearby van Hove sin- 
gularities (vHS). The physics in this case becomes more 
transparent if the total current is split into the contribu- 
tions due to the different subbands. This is achieved by 
resorting to a unitary transformation which diagonalizes 
the electronic Hamiltonian for each layer perpendicular 
to the transport direction (s^j- For a (iV, 0) nanotube, 
the different subbands correspond to linear chains with 
alternating hoppings 70 and 7^ = 270 cos((j'7r/A^) {q = 
0,1...,A^ — 1). When N is an integer multiple of three, 
the metallic subbands correspond to q = N/3, 2N/3. Fig- 
ure 2 shows the current pumped as a function of the gate 
voltage Vgate (inset) and driving frequency (main frame) 
applied to a (24,0) nanotube for a typical set of param- 
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FIG. 2: Maximum current pumping obtained for (24,0) car- 
bon nanotube. (color online) Inset: Pumped current as a 
function of the gate voltage. The different curves are the con- 
tributions from the different sub-bands. The position of the 
van Hove singularities are indicated with arrows on the hori- 
zontal axis. The driving frequency used in the inset is marked 
with a vertical dashed line in the main frame. Main frame: 
Maximum pumped current for each contributing sub-band as 
a function of the driving frequency. el4c is set as in Fig. 1 
and L — 100 nm. 



eters. Different curves show the contributions coming 
from each subband with q = 7,8,9 to the total current 
(since the contributions for q ~ 15, 16, 17 are equal to 
those shown here we concentrate only in the latter) . The 
curve labeled as g = 8 corresponds to one of the two lin- 
ear bands which crosses the CNP and which is responsi- 
ble for the metallic character of the tube. The remaining 
curves, q = 7 and q ~ 9, lead to the first and second 
van Hove singularities, respectively. By inspection of the 
inset one sees that subbands labeled with q = 7,9, which 
reveal vHSs within the displayed gating range, dominate 
over the one with q = 8. Furthermore, the contribution 
due to q = 7,9 show a damping behaviour, revealing a 
maximum intensity close to the corresponding vHS and 
then a decay. What triggers such a dramatic increase 
of the current shown in Fig. 2 inset? Careful analysis 
of the numerical results show that as the gate voltage 
is tuned closer to the vHS, there is a crossover from a 
regime where inelastic effects (which are responsible for 
non-adiabatic pumping) are weak to one where they be- 
come dominant (close to the vHS where the level spacing 
becomes smaller than the driving frequency hfl ). 



Another relevant point is the scaling of the pumped 
current with the driving frequency. In contrast to 
adiabatic pumping, where the frequency can be made 
very small while keeping a constant pumped charge 
per cycle, for non-adiabatic pumping through a system 
with well separated resonances, the frequency scaling of 
the pumped current is typically found to be quadratic 
[1: El ^M- the main frame of Fig. 2, we show the 
frequency scaling of the maximum pumped current (dis- 
criminated by subband) in a broad range of frequencies 
going from tens of GHz to THz. Although, no simple 
behavior is expected for the maximum current due to a 
complicated interplay between finite size effects and driv- 
ing parameters. Fig. 2 shows a dependence on frequency 
for q = 7,9 which is for lower frequencies definitely bet- 
ter than the quadratic one observed for q — 8. Notwith- 
standing, all the curves show a decrease for high enough 
frequency. This is shown only for g = 9 in the plot where 
the current reveals a maximum value. This behavior can 
be understood as an interplay between the enhancement 
of inelastic processes for low frequency as the vHS be- 
comes prominent and the opposite tendency produced 
when the frequency is large enough such that the field be- 
comes ineffective in exciting photons. The same pumping 
behaviour was verified for armchair graphene nanorib- 
bons, the only difference being the precise position of the 
vHS and the number of available channels. 

Conclusions. Most of the studies on quantum charge 
pumping focus on the behavior of a two-parameter pump 
close to isolated resonances. Here wc focused on a differ- 
ent scenario where non-adiabatic effect dominate: single- 
parameter pumping in a system with very good contacts. 
Our results are illustrated for two cases of much ex- 
perimental relevance: carbon nanotubes and graphene 
nanoribbons. Besides offering currents that are higher 
than those driven close to the charge neutrality point, 
pumping phenomena held close to a van Hove singularity 
renders an improved scaling with frequency. On the other 
hand, the pumped current close to the charge neutrality 
point may also be substantially enhanced by tuning the 
driving parameters and system size. 
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